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ABSTRACT 

We present the analysis of MOA-2007-BLG-197Lb, the first brown dwarf companion to a Sun-like star detected through gravitational 
microlensing. The event was alerted and followed-up photometrically by a network of telescopes from the PLANET, MOA, and 
pFUN collaborations, and observed at high angular resolution using the NaCo instrument at the VLT. From the modelling of the 
microlensing light curve, we derived basic parameters such as, the binary lens separation in Einstein radius units (s ^ 1.13), the mass 
ratio q = (4.732 ± 0.020) x 10 -2 and the Einstein radius crossing time (f E == 82 d). Because of this long time scale, we took annual 
parallax and orbital motion of the lens in the models into account, as well as finite source effects that were clearly detected during the 
source caustic exit. To recover the lens system’s physical parameters, we combined the resulting light curve best-fit parameters with 
( J,H,K S ) magnitudes obtained with VLT NaCo and calibrated using IRSF and 2MASS data. From this analysis, we derived a lens 
total mass of 0.86 ± 0.04M 0 and a lens distance of D L = 4.2 ± 0.3 kpc. We find that the companion of MOA-2007-BLG-197L is a 
brown dwarf of 41 ± 2Mj observed at a projected separation of a ± = 4.3 ±0.1 AU, and orbits a 0.82 ± 0.04 M 0 G-K dwarf star. We 
then placed the companion of MOA-2007-BLG-197L in a mass-period diagram consisting of all brown dwarf companions detected 
so far through different techniques, including microlensing, transit, radial velocity, and direct imaging (most of these objects orbit 
solar-type stars). To study the statistical properties of this population, we performed a two-dimensional, non-parametric probability 
density distribution fit to the data, which draws a structured brown dwarf landscape. We confirm the existence of a region that is 
strongly depleted in objects at short periods and intermediate masses (P < 30d, M ~ 30 - 60 Mj), but also find an accumulation of 
objects around P ~ 500 d and M ~ 20 Mj, as well as another depletion region at long orbital periods ( P > 500 d) and high masses 
(M > 50 Mj). While these data provide important clues on the different physical mechanisms of formation (or destruction) that shape 
the brown dwarf desert, more data are needed to establish their relative importance, in particular as a function of host star mass. 
Future microlensing surveys should soon provide more detections, in particular for red dwarf hosts, thus uniquely complementing the 
solar-type host sample. 
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1. Introduction 


Gravitational microlensing is a powerful technique for detecting 
extrasolar planets ( jMao & Paczynski||1991[ >, and it holds great 
promise for detecting populations of brown dwarf companions 
to stars. Compared to other detection techniques, microlensing 
provides unique information on the population of exoplanets, be¬ 
cause it allows the detection of very low-mass planets (down to 
the mass of the Earth) at long orbital distances from their host 
stars (typically 0.5 to 10 AU). It is also the only technique that 
allows discovery of exoplanets and brown dwarfs at distances 
from the Earth greater than a few kiloparsecs, up to the Galactic 
bulge, which would have been hard to detect with other methods. 

Exoplanets are found to be frequent by all detection tech¬ 
niques (e.g., Cassan et al.|2012[|Bonfils et al.|2013)|Mayor et al.| 


|2011||Sumi et al.|2011||Gouldetal.|2010b| ), and recent statistical 

microlensing studies even imply that there are, on average, one 
or more bound planets per Milky Way star (jCassan et al.|2012|>. 


Conversely, brown dwarfs appear to be intrinsically rare, to the 
point that shortly after the first exoplanet detections, it led to the 
idea of a “brown dwarf desert” ( |Marcy & Butler||2000) > bridg¬ 
ing the two well-defined regions of binary stars and planetary 
systems. While in the past, brown dwarfs were defined as ob¬ 
jects of mass within the deuterium- and hydrogen-burning limits 
(13 - 74 Mj, Burrows e t al|200T) , it appears today that differ¬ 
ent formation scenarios can build objects with similar masses 
but with different natures (super-massive planets, or low-mass 
brown dwarfs). An object formed via core accretion and reach¬ 
ing 13 Mj would, for example, be able to start deuterium burn¬ 
ing, as would an object of same mass formed by gravitational 
collapse of a cloud or in a protoplanetary disk ( Molliere & Mor-[ 
|dasini l2Ql2] >. 


Despite their low occurrence, a number of brown dwarf com¬ 
panions to stars have been discovered by different methods: ra¬ 
dial velocity and transit (e.g., Moutou et al. 2013, Dfaz et al. 


120131 [Sahlmann et al.||20lTt|Johnson et al.|201 1| |Deleuil et al 


Article number, page 1 of 


16 

































































































































A&A proofs: manuscript no. draft 


2008 ), direct imaging (e.g., |Lafreniere et al.| |2007| ) and mi- 
crolensing. With regards to microlensing, there are still a few 
brown dwarf detections, but this is mainly because until now 
observing priority has been given to exoplanets. Nevertheless, 
these detections are of particular interest because they provide 
a unique view of brown dwarfs around low-mass stars (mainly 
M dwarfs), which complements the currently available sample 
mostly composed of solar-type stars. 

Microlensing detections so far include isolated brown 
dwarfs, brown dwarfs hosting planets, and brown dwarf com¬ 
panions to stars. The first isolated brown dwarf detected through 
gravitational microlensing (a 59 ± 4Mj brown dwarf located 
at 525 ± 40 pc in the thick disk of the Milky Way) was re¬ 
ported by Gould et al. 1(2009) in microlensing event OGLE-2007- 


BLG-224. Two brown dwarfs with planetary-mass companions 
were discovered in events OGLE-2009-BLG-151/MOA-2009- 
BLG-232 and OGLE-201 l-BLG-0420 by |Choi et al. ( 2013) . In 
both cases, the planets were super Jupiters (7.9 ± 0.3 Mj and 
9.9 ± 0.5 Mj, respectively) with the hosts being low-mass brown 
dwarfs (19 ± 1 Mj and 26+1 Mj, respectively), with very tight 
orbits (below 0.4 AU). Similarly, |Han et aF ( 2013[ ) report a 
23 ± 2 Mj field brown dwarf hosting a 1.9 + 0.2 Mj planet in 
a tight system after the analysis of the event OGLE-2012-BLG- 
0358. 

The first published microlensing detection of a brown dwarf 
companion to a star is OGLE-2008-BLG-510/MOA-2008-BLG- 
369, which was first reported by Bozza et al. ( 2012[ > as an am¬ 
biguous case between a binary-lens and a binary-source event. 
The data were reanalysed by Shin et ah] ( |2012a| ), who con¬ 
cluded that the binary-lens model involving a massive brown 
dwarf orbiting an M dwarf was preferred. |Shin et akj ( 2012b) 
conducted a database search for brown dwarf companions by 
focusing on microlensing events that exhibit low mass ratios. 
Among seven good candidates with well-determined masses 
(combination of Einstein radius and parallax measurements), 
they found two events that involve brown dwarfs: OGLE-2011- 
BLG-0172/MOA-201 l-BLG-104 with mass 21 ± 10 Mj around 
an M dwarf, and MOA-201 l-BLG-149, a 20 + 2 Mj brown dwarf 


also orbiting an M dwarf. Similarly, Bachelet et al.j (2012 1 re¬ 
ported the detection of another ~ 52 Mj brown dwarf orbiting 
an M dwarf in MOA-2009-BLG-411, although the lens mass 
could not be determined exactly, and was estimated through sta¬ 
tistical realisations of Galactic models. In microlensing event 
MOA-2010-BLG-073, [Street et al.] ( |20T3] l find that the lens 
was composed of a 11.0 ± 2.0 Mj companion (hence near the 
planet/brown boundary) orbiting an M dwarf of 0.16 ± 0.03 M 0 . 


Jung et al. (2015) report the detection of a star at the limit of the 


brown dwarf regime hosting a companion at the planet/brown 
dwarf boundary (13 ± 2Mj). More recently, [Park et al ] ( |2015[ > 
have reported the discovery of a binary system composed of a 
33.5 ± 4.2 Mj brown dwarf orbiting a late-type M dwarf in mi¬ 
crolensing event OGLE-2013-BLG-0578. 

Here we report the first microlensing discovery of a brown 
dwarf orbiting a Sun-like star. This new brown dwarf has a mass 
of 42 Mj and it was observed at a projected separation of 4.3 AU 
from its G-K dwarf host star. 

In sec. 2, we present photometric data collected on MOA- 
2007-BLG-197L by several round-the-world telescopes, as well 
as high resolution adaptative optics images taken with NaCo at 
VLT. In sec. 3, we present the full analysis of the light curve and 
in sec. 4 we derive the physical parameters of the lens by com¬ 
bining all independent informations. In sec. 5, we adopt a statis¬ 
tical point of view to analyse the current population of brown 
dwarfs detected with different methods using non-parametric 
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probability density estimation tools. In sec. 6 we summarise our 
results, and underline the importance of future microlensing ob¬ 
servations to characterise the populations of objects in the mass 
region between planets and stars. 


2. Observational data 

2.1. Alert and follow-up 

MOA-2007-BLG-197 (/ = 359.711, b = -5.509, or RA (J2000) 
= 18:07:04.729, DEC (J2000) = -31:56:46.77) is a microlens¬ 
ing event that was alerted by the MOA collaboration (1.8m tele¬ 
scope located in Mount John, New Zealand) in 2007 May 28 (or 
THJD - 424 l Q on the photometric light curve shown Fig. [TJi. 
Soon after, between 1-5 June (THJD - 4251-4255), a very 
cloudy weather in New Zealand seriously affected the quality 
and reliability of the photometry, resulting in gaps in the time 
series. A magnification peak in the MOA light curve was passed 
around June 6 (THJD =* 4258). On June 8 (THJD =* 4259.7), 
the PLANET collaboration added the event to its target list as 
a regular mid-magnification object. Follow-up observations in¬ 
clude data from the PLANET network using the Danish tele¬ 
scope 1.54m at La Silla (Chile), Canopus lm in Hobart (Tas¬ 
mania), SAAO lm in Sutherland (South Africa) and Perth 0.6m 
telescope (Australia), while the //FUN collaboration collected 
data from CTIO 1.3m in Mount Cerro Tololo (Chile). 

On June 19 (THJD - 4272), PLANET observers noticed 
that a few MOA data points were slightly above the standard 
single-lens theoretical curve. As this usually happens quite often 
in microlensing data, no alert was released. Indeed SAAO soon 
revealed that the alleged deviation was a false alert, but follow¬ 
up observations were continued. On July 5 (THJD ^ 4287), 
even though the full moon affected the quality of the observa¬ 
tions, a public alert was issued after Danish 1.54m data were 
found to be above the single lens curve by more than 0.2 mag 
for more than five consecutive days. Intensive follow-up obser¬ 
vations from SAAO confirmed the rise in brightness, announcing 
a caustic crossing. While Perth was overclouded, Canopus then 
took over, and was the only telescope to densely cover the caus¬ 
tic exit, which took place during the night of July 4 in Australia 
(THJD =/ 4287.0-4287.4, see inset of Fig. |T}. 

2.2. Photometric light curve 

The photometric data was reduced several times using different 
software to check their consistency. The final lightcurve data 
from all telescopes were extracted with the PLANET pipeline 
PySIS ( |Albrow et al.|[2009] >, which is a DIA-based algorithm 
(Difference Image Analysis, |Alard|2000||Alard & Lupton|l998) 1 . 
For images taken in particularly bad weather conditions, we ex¬ 
amined by eye each image in order to check whether the subtrac¬ 
tion was correct or not. In the latter case we had to exclude them, 
but on the basis of image quality only. In the post-processing 
of the reduced data, we applied a cut in seeing and sky back¬ 
ground, although in a very conservative way so that possible 
low-amplitude signals were not rejected. Not surprisingly, most 
of the data taken by MOA during the very bad weather period 
mentioned before (THJD - 4251-4255) had to be discarded. 
Unfortunately, a critical light curve feature is thereby not cov¬ 
ered (which has important consequences, see Sect. [3]), however 
there was no other choice but to remove these data, as they might 
otherwise have affected the reliability of our models. 

1 THJD = HJD - 2,450,000. 
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Table 1. Telescopes and photometric data sets. 


Telescope 

Location 

Filter 

Data 0 

f b 

MOA (1.8m) 

Mount John, New Zealand 

Rm c 

504 

1.4 

PLANET Danish (1.54 m) 

La Silla, Chile 

I 

167 

1.7 

PLANET Canopus (1.0 m) 

Mount Canopus, Tasmania 

I 

45 

3.0 

PLANET SAAO(l.Om) 

Sutherland, South Africa 

I 

43 

1.6 

yuFUN CTIO (1.3 m) 

Mount Cerro Tololo, Chile 

I 

28 

1.3 

PLANET Perth (0.6 m) 

Perth, Australia 

I 

15 

1.4 


Notes. Number after data cleaning. (i) Error bar rescaling factor (Sect. 2.2 '. (c) MOA broad R/I filter. 



4220 4240 4260 4280 4300 4320 4340 4360 

HJD- 2,450,000 


Fig. 1. In the upper panel, the light curve of MOA-2007-BLG-197 and the best-fit model (solid line) are plotted with a zoom on the caustic exit 
on the right-hand side. On the left-hand side, the structure of the resonant caustic is drawn in red, as well as the trajectory of the source, in black 
(axes are in Einstein radius units). The source is too small to be distinguished. The four intervals £i _4 indicate time intervals where possible source 
caustic crossing (caustic entry) have been investigated. In the lower panel, the residuals of the best-fit model are shown. 


The final light curve data amounts to a total of 802 data 
points. They are summerised in Table |T] As seen in the table, 
all telescopes use a similar I -band filter, apart from the MOA 
1.8m telescope which is equipped with a broad R/I filter (re¬ 
ferred to as Rm)- Additionally, a few V-band images were taken 
by PLANET and /iFUN to produce colour-magnitude diagrams 
(see Sect. |4.1| >. 

The last concern about the photometry was the estimation 
of the error bars of the data. Galactic bulge fields are highly 
crowded with stars, and during a microlensing event, the flux 
variation can easily span two order of magnitudes for high- 
magnification events. These pose severe challenges to manage 
a good estimation of the error bars. As a matter of fact, in mi¬ 
crolensing experiments it is long known that data reduction soft¬ 
ware usually underestimates error bars. Furthermore, error bars 
can vary significantly from one data set to another, with the risk 


that one data set dominates over the others at the modelling 
stage. A relatively robust method to prevent these drawbacks is 
to rescale the error bars, based on the best model fitting the data. 
For each data set, the (classical) x 1 is set up to the number of de¬ 
grees of freedom by adjusting a rescaling factor / in the formula 
cP 2 = f~cr~ + cr 2 , where cr' and cr are respectively the rescaled 
and initial error bars on the magnitudes, and <x 0 = 4 x 10 4 a 
constant accounting for the data most highly magnified. The / 
factors are given in Table [1] for each data set. 

2.3. VLT NaCo high resolution images 

On the night 20/21 of August 2007 (THJD =* 4333.0) we ob¬ 
tained first epoch observation^] of high resolution adaptive op- 

2 ESO Programme ID 279.C-5044(A). 
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tics (AO) images in the near-infrared bands J, H and K s using 
the NaCo instrument, mounted on 8.2m ESO VLT Yepun tele¬ 
scope (Fig. [2). The source star was then still magnified by a fac¬ 
tor of about A - 1.30. A year later, on the night 2008 August 
3/4 Q’HJD - 4682.1), we carried out second epoch observa- 
tionsnwhile the event was back to baseline magnitude. In princi¬ 
ple, when two epochs are obtained at different magnifications, it 
should be possible to directly disentangle the lens flux from the 
source flux. In our case, however, the combination of a relatively 
high blending factor and low magnification did not support this 
direct measurement. Therefore, in the global analysis we used 
only the first epoch images. We reduced and calibrated the NaCo 
images following the general method outlined in [Kubas et al.[ 
(20121 and briefly described below. 



NACO 


IRSF * 


Fig. 2. The image on the left shows a sub-region (94" x 77") of an 
original 3' x 3' A - ,-band IRSF image used to calibrate the NaCo magni¬ 
tudes. The image on the right is the corresponding A'j-band NaCo image 
(27" x 25") used to cross-identify the stars. 


magnitude^jof star 9160-S870-T in the listed standards of Pers- 
son et al. ( 1998) , which had the following advantages: this star 
was brighter than the previous 2MASS reference, and the smaller 
pixel scale made it less sensitive to blending contamination in 
our crowded field. For consistency, however, we measured the 
magnitude of 2MASS 18070520-3156409 in the NaCo frame, 
and found an agreement to better than 3%. 

The second step was to extract accurate photometry from AO 
images in the infrared. This was not a straightforward task, since 
the shape of the point-spread function (PSF) often is not well 
fitted by analytical profiles, and also depends on the position of 
the target with respect to the star used for AO correction. Follow¬ 
ing ^Sbas£Tar](j20l2|), we constructed a PSF reference directly 
from stars in the NaCo frame, using the StarFinder package 
olaiti et al. 2000). This software was especially designed for AO 


images of crowded stellar fields. We found J = 17.68 ± 0.06, 
H = 17.05 ± 0.05 and K s = 16.87 ± 0.05. The quoted error 
bars are dominated by the uncertainties in determining the true 
PSF shape and the scatter of the sky and unresolved background 
sources. The measurements obtained with StarFinder are sum¬ 
marised in Table [2] 


Band 

Epoch 

Magnitude 

Date [THID] 

FWHM 

J 

1 

17.68 + 0.06 

4333.03906250 

0.19” 

J 

2 

17.67 ± 0.05 

4682.12500000 

0.17” 

H 

1 

17.05 ± 0.05 

4333.05468750 

0.14” 

H 

2 

17.04 ± 0.04 

4682.14453125 

0.11” 

K s 

1 

16.87 + 0.05 

4333.02343750 

0.11” 

K s 

2 

16.89 ± 0.04 

4682.10546875 

0.11” 


Table 2. Derived NaCo magnitudes of the microlensing event target in 
the (J, H, K s ) filters for the two epochs 2007 August 20/21 (epoch 1) and 
2008 August 3/4 (epoch 2). The values of Epoch 1 are used to constrain 
the lens mass-distance relation. 


The night was clear and stable according to the observatory 
night logs, and the target was observed at low airmass and in 
good seeing conditions (< 0.8")- The data were taken in auto¬ 
jitter mode within a 10" jitter box including the target in order 
to be able to correct for bad pixels and sky background. The 
AO correction was done using a K - 11.17 ± 0.02 mag star 
(2MASS 18070464-3156423) at an angular separation of about 
5" from the microlensing target. The data were dark subtracted, 
flat-fielded and co-added with the tools underlying the NaCo 
pipeline software |Devillard|l999) . 

To derive the photometry from the reduced data, the first step 
was to compute the zeropoints for the conversion of instrumental 
magnitudes to calibrated magnitudes. A first possibility was to 
use catalogued stars within the target frame field-of-view (FOV), 
in our case 28" x 28". While several cross-matches between the 
NaCo frames and the 2MASS catalogue were identified, only 
2MASS 18070520-3156409 (J = 14.24+0.05, H = 13.69+0.05, 
K = 13.48 ± 0.05) turned out as suitable (other potential calibra¬ 
tors in the NaCo images were either saturated or in the nonlinear 
regime of the detector). We finally checked that this star was not 
variable, by comparing a series of //-band images taken with 
Andicam at CTIO, which is well calibrated to 2MASS thanks 
to its 2.4' x 2.4' FOV. A second option for calibration was to 
use the zeropoints derived from the photometric standards taken 
with NaCo directly before the observations on the night 20/21 of 
August and at similar airmass. For calibration we used J, H , K 


ESO Programme ID 381.C-0425(A). 
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The final step consisted in correcting the target for inter¬ 
stellar extinction, by fitting the position of the red clump giants 
(RCG) in the three colour-magnitude diagrams (CMD) involving 
the measured /, H, K s reddened magnitudes. Because the nomi¬ 
nal range in magnitudes from NaCo is above the 2MASS faint 
limit, we performed InfraRed Survey Facility (IRSF) observa¬ 
tions at SAAO to extend the available 2MASS star list into the 
regime of stars measured within the NaCo frame. We used IKato] 
|et al.|(|2007) to obtain the calibration of IRSF images with re¬ 


spect to the 2MASS reference-star catalogue, noting that Janczak 
|et al.| ( |20T0) found that no additional colour term is needed be¬ 
tween NaCo and IRSF filters (the full process to build the cal¬ 
ibration ladder is detailed in jKubas et al.||20l 2) . Data to con¬ 
struct the NaCo+IRSF CMD are extracted from 2304 stars iden¬ 
tified within a 3' circle around the target in IRSF images, and 
135 stars identified in the NaCo images (Fig. [2). The result¬ 
ing de-reddened and calibrated CMD is plotted in Fig. [3] with 
a fit of the Red Clump Giant (RCG) position using a 10 Gyr, 
Z = 0.019 isochrone from Bressan et al. (2012), and assuming 
a RCG distance modulus of p — 14.6 mag. This distance mod¬ 
ulus is longer than the value used in |Kubas et al. (2012) since 
the event is located in the Galactic bar (see Sect. 4.1 1 for a de¬ 
tailed discussion). From this analysis, we derived the following 
extinctions: Aj = 0.51 ± 0.05 mag. Ah = 0.33 ± 0.05 mag and 
Ak s = 0.22 ± 0.05 mag. These data are used in Sect.|4.2|to con¬ 
strain the lens and mass distance. 


4 Note that the difference in the transmission profile between K and Ks 
bands is less than 1%, so negligible in the present case. 
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2.0 


Fig. 3. (i, J—K s ) colour-magnitude diagram in the 2MASS system com¬ 
bining NaCo and IRSF data (respectively, in blue and black). The red 
dot is the magnitude/colour of the source derived from the source char¬ 
acterisation (Sect. |4~T| i. The values are not corrected for interstellar ex¬ 
tinction. The red clump (RCG) is fitted by the red over-plotted 10 Gyr, 
solar-metallicity isochrone from |Bressan et al7] ( |201 2) . 

3. Light curve modelling 

3. 1. Binary lens parameters 

We start by modelling the light curve of MOA-2007-BLG-197 
with a static binary-lens model, in order to identify broad classes 
of possible solutions. In this model, the lens is characterised by 
its binary mass ratio q, and ,v, the projected separation of the 
binary lens in angular Einstein radius units 0 E . Here, 


6 e - 


AGMlDs \ 
c 2 D s \D L )’ 


(1) 


where G is the gravitational constant, c the speed of light, De 
and Ds are respectively the observer-lens and observer-source 
distances and M the total mass of the lens. Four additional pa¬ 
rameters describe the source mean rectilinear motion: the mini¬ 
mum impact distance uq between the source and the origin (here, 
the centre of mass of the system, with the more massive body on 
the right-hand side), to, the date at which the source reaches mq, 
f E , the time it takes for the source to travel one Einstein angular 
radius, and a, the angle between the source trajectory and the 
lens symmetry axis. 

When the source approaches a caustic, finite-source effects 
cannot be neglected and substantial deviations from a point- 
source model are expected. Finite-source effects are included in 
the modelling through 


P = 


Os 

0e 


( 2 ) 


where 6s is the angular radius of the source, and p the same quan¬ 
tity but in Einstein radius units. Here and in the following, a 
point-source model is used when the source is far enough from 
the caustics (t < 4240.0 and t > 4300.0); closer to the caustics, 


a hexadecapole approximation to the finite-source magnification 
( |Gouldl2008[ ) is used (t e [4240.0,4247.9] U [4260.0,4286.0] U 


[4287.5,4300.0]), until it breaks down, approximately when the 
source is closer than ~ 3 p from the caustics. In this case, a full 
integration along the source images contours (Dominik 2007; 
Bozza||2010[ Gould & Gaucherel 1997 1 or ray-shooting of the 


images are required (Dong et al. 2006| |2009| |Bennett||2010| > 
to compute the magnification. Since finite-source effects dra¬ 
matically increase the computational cost, these time intervals 
are reduced as much as possible (here, t e [4247.9,4260.0] U 
[4286.0,4287.5]). 

Microlensing light curves are usually sensitive to stellar limb 
darkening (e.g., |Albrow et al.||1999b) |Cassan et al. |2006|). We 
thus model the source as a linear limb-darkened disk ( |An et al.| 
2002 [ Zub et al. 201 1]> described by the intensity-normalised pro¬ 


file 
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i -r 



(3) 


where 0 < r < 1 is the fractional radius, and T the linear 
limb-darkening (LLD) coefficient. While in special cases limb- 
darkening laws beyond the linear law may slightly improve the 
model (e.g., |Cassan et al.|2004[|Kubas et al.|2005) , in the case of 
MOA-2007-BLG-197L (relatively sparse data coverage of the 
caustic crossings) it is an excellent approximation. As seen in 
the right inset of Fig. |T] only Canopus data are sensitive to limb 
darkening, although the light curve sampling is not dense enough 
to provide strong constraints on the LLD coefficient. We there¬ 
fore use Claret & Bloemen ( |201 1] > for the source surface grav¬ 
ity and effective temperature found in Sect. 14. 1 1 (log g ~ 4.47, 
r eff ~ 5350 K) to estimate MOA-2007-BLG-197L LLD coeffi¬ 
cients, and found ~ 0.4 for the / filter and ~ 0.5 for the R filter. 
We checked that refining the Canopus-/ LLD coefficient with a 
fit to the data leads to L/ = 0.48. We then adopted E/ = 0.48 for 
all /-band data and Tr m = 0.5 (MOA Tr m filter is a broad R/I 
filter). 

Finally, two additional parameters, F' s (source flux) and F' H 
(blended flux), are included per individual observatory or filter, 
so that for a given individual data set i the total flux of the mi¬ 
crolensing target reads 


F‘(t) = A(t) F' s + F l R , 


(4) 


where A(t) is the time-dependent source flux magnification fac¬ 
tor. The blending flux Cg accounts for all luminous contribu¬ 
tions other than the source, in particular, it includes the flux from 
the lens. For the best model presented in Table [3] the blending 
flux ratios L’b/T's are respectively 17.79 (MOA), 38.14 (SAAO), 
12.85 (CTIO), 3.97 (Canopus), 66.97 (Perth), 12.10 (Danish). 
The high value for Perth is easily explained by the fact that the 
aperture of the telescope is small and the seeing was relatively 
high, which results in a higher blend fraction. 


3.2. Caustic-crossing parametrisation 

The shape of the MOA-2007-BLG-197 light curve (Fig. [TJ 
clearly indicates that the source crosses (i.e., exits) a fold caus¬ 
tic at THID =* 4287.2, mainly thanks to Canopus data densely 
covering this part of the light curve. On the other hand, the date 
of the caustic entry is ambiguous. In fact, for such a small angu¬ 
lar source size as suggested by the caustic exit, the caustic entry 
could easily fit in one of the several gaps in the data coverage, in 
particular those marked as £ 1-4 in Fig. [T] 
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One of the modelling challenges for this event was to find 
all possible local minima involving all possible caustic entry 
dates. The problem with the set of parameters (mq, a, t E , to) de¬ 
scribed in the previous section is that it is not well suited to 
explore efficiently the parameter space. In fact, only very spe¬ 
cial combinations of these parameters produce caustic crossings 
at the right locations on the light curve. |Albrow et al.| ( |l999a[ l 
proposed for the first time to introduce specific parameters to 
model a caustic crossing. Cassan ( 2008j > developed further this 
approach, and generalised it to a pair of caustic crossings (entry 
and exit) by introducing a specific parametrisation of the caus¬ 
tic curves. This method is particularly efficient for an event like 
MOA-2007-BLG-197. 

In this formalism, two alternative parameters and f out are 
used to fit for the dates of the source entry (in) and exit (out), 
while two other parameters ,v ln and s out (which are curvilin¬ 
ear distances along the caustic curve) are used to fit for the 
source centre ingress and egress points on the caustic. It has 
been demonstrated that this parametrisation is more efficient in 
locating all possible fitting source-lens trajectories ( |Kains et al.| 
2009 2012), because they all produce caustic-crossing features 


at the observed dates. Cassan et al. j( |2010| l later derived Bayesian 
priors on parameters (f; n , f out , ,v in , s out ) to explore even more effi¬ 
ciently the parameter space. The best-fit values derived from the 
minimisation process are finally converted back to the classical 
parameters. 

We explorec0 all static binary-lens models using this 
parametrisation for all possible source caustic crossing dates 
within the £j to £4 intervals displayed in Fig.[T| and for a regular 
(log s, log q) grid of 30 x 30 spanning log(0.4) < log s < log(4) 
and -4 < log q < 0. The minimisation was performed us¬ 
ing a classical Markov Chain Monte Carlo (MCMC) algorithm. 
The best model involves a caustic entry in the time interval £ 1 , 
the longest of the four gaps explored. The best-fit values of the 
model parameters and the corresponding y 2 are listed in column 
ESBL (Extended-Source, Binary-Lens) of Table[3] and posterior 
probability densities (correlation plots) are presented in Fig. [9] 


3.3. Annual parallax 

Given the long timescale found previously for the static binary 
lens model (f E ~ 80 days), this event is a priori likely to exhibit 
annual parallax effects. The relative lens-source (annual) paral¬ 
lax is expressed as 


Tfiel - 


AU 


AU 
A ’ 


so that the Einstein angular radius Eq. <jT]) also reads 


9 e = ^kM 7i , 


(5) 


( 6 ) 


where k - 8.144 mas/M Q . The Einstein parallax vector 7 t/. has 
an amplitude 


7T E = 


^l-el 


(7) 


and a direction along the lens-source proper motion. As fitting 
parameters, it is convenient to decompose tie into a northwards 
component 7Te,n and an eastwards component 7Te,e (|An et al. 
|2002| l, with advantages discussed in |Gould| ( |2004| l. 


5 These very intensive computations were performed on the Tasmanian 
cluster TPAC (Tasmanian Partnership for Advanced Computing). 
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The best-fit model including parallax only (besides the pa¬ 
rameters described in Sect. ED is presented in column ESBL+P 
of Table [3j and correlation plots are shown in Fig. [10] It can be 
noticed that the binary mass ratio q and separation s change lit¬ 
tle by including parallax compared to ESBL, resulting in a very 
similar resonant caustic structure. As expected, including par¬ 
allax in the model improves the y 2 , by around 100. Neverthe¬ 
less, we found that several solutions with very different values 
of parallax (between 0.8 to 1.5) give almost identical y 2 differing 
by only one or two units. The differences between these models 
come from differences in the caustic entry date t m , i.e. precisely 
where no data are available to constrain the model. Hence, al¬ 
though parallax improves the fit, it is unlikely that a workable 
measurement can be obtained (we discuss this further later). 


3.4. Orbital motion 

Orbital motion of the lens is also a priori likely to produce no¬ 
ticeable effects, because the event is fairly long and the caustic 
resonant. The orbital rotation of the two lens components affects 
the caustic in two ways: it changes the projected binary separa¬ 
tion s and the orientation of the caustic in the plane of the sky. To 
first order, one can write s(t) - s+s {t—t r ) and a(t) - a+a (t-t r ), 
where f, is an arbitrary reference date chosen close or equal to to. 
These two effects are included in the model through parameters 
yy - s/s and y ± = a, with y 1 = + y]_. 

The best-fit model including lens orbital motion with ESBL 
is presented in column ESBL+LOM of Table[3](correlation plots 
are shown in Fig. m According to this model, the changes in 
the caustic geometry are very slow, which should result in a poor 
signature on the light curve. Therefore, while this model also 
provides a better fit than ESBL alone (Ay 2 = 93), it is unlikely 
to give a good constraint on the lens orbital motion. Further¬ 
more, the Ay 2 between this model and ESBL+P (parallax alone) 
is only lower than 10 although two more parameters have been 
included in the model. We conclude that the key time intervals 
for disentangling orbital motion from parallax are not covered 
densely enough by the data. 


3.5. Annual parallax and orbital motion 

As (annual) parallax and orbital motion individually lead to a 
similar improvement of the ESBLy 2 , we now include both ef¬ 
fects in the model to check whether a bettery 2 can be found. Par¬ 
allax and orbital motion are known to be very correlated (Batista 
|et al.| 20 fT]|Skowron et al.| 2011 j i. 

We investigate the two cases Uq > 0 and uq < 0 (columns 
ESBL+P+LOM of Tabled to check for possible ecliptic degen¬ 
eracy (Skowron et al. 2011[ >, and find that uq > 0 is preferred 
by Ay- = 11.9. The ESBL best-fit parameters, again, are rela¬ 
tively stable. The correlation plots of the best model are shown in 
Fig. [12] As expected, the overall fit is not significantly better than 
ESBL+P (Ay 2 = 12.2) or ESBL+LOM alone (Ay 2 = 20.5). This 
is consistent with a strong degeneracy between these parameters, 
that cannot be reliably broken in the case of MOA-2007-BLG- 
197. 


3.6. Summary and conclusions from the modelling 

A model with (annual) parallax improves the goodness-of-fit by 
Ay 2 ~ 100 compared to a static binary lens, but several models 
with very different values of 7r E give comparable values ofy 2 . 
Moreover, models with parallax or orbital motion alone lead to 
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Table 3. Best-fit solutions for the different models of MOA-2007-BLG-197. 


Parameter [unit] 



Model 



ESBL n 

ESBL + P* 

ESBL + LOM c 

ESBL + P + LOM 

mo < 0 uo > 0 


;r/d.o.f. 


Ar =x 2 -xLt 

S 

ql 10- 2 

t E [days] 

Mo/10 2 
to [THJD] 
a /10 1 [rad] 
p/io- 4 

TtE.N 

7tE,E 

s/lO -1 [rad/year] 
d/lO -1 [rad/year] 


(942.2/795) = 1.19 

113.5 

1.1291 ±0.0012 
4.972 ± 0.043 
80.72 ±0.21 
6.39 ±0.10 
4259.31 ±0.10 
9.123 ±0.056 
5.69 ± 0.23 


(840.9/792) = 1.06 
12.2 

1.1347 ±0.0020 
4.867 ± 0.080 
78.39 ±0.61 
5.57 ±0.15 
4259.72 ± 0.40 
8.54 ± 0.20 
5.91 ±0.30 
0.82 ± 0.25 
-0.31 ±0.15 


(849.2/793) = 1.07 

20.5 

1.1620 ±0.0044 
3.48 ± 0.25 
79.4 ± 1.0 
7.17 ±0.32 
4258.34 ± 0.70 
7.57 ±0.17 
5.20 ± 0.27 


1.22 ±0.16 
0.24 ± 0.65 


(840.6/790) = 1.06 
11.9 

1.1361 ±0.0086 
4.934 ± 0.044 

82.6 ± 1.4 
-5.43 ±0.12 
4259.23 ±0.10 
-9.62 ±0.14 
6.13 ±0.33 
1.32 ±0.31 
0.82 ± 0.65 
-0.04 ±0.16 
-2.76 ± 0.60 


(828.7/790) = 1.05 

0 

1.1254 ±0.0044 
4.732 ± 0.020 
82.3 ± 1.2 
5.59 ± 0.28 
4258.884 ±0.013 
9.709 ± 0.060 
5.30 ± 0.24 
-0.45 ±0.13 
-0.96 ±0.15 
-0.46 ±0.13 
2.38 ± 0.60 


Notes. <fl) Extended-source binary-lens model. <A) Microlensing parallax. <c) Lens orbital motion. 


comparable^ 2 . This is not surprising for two reasons. Firstly, the 
caustic entry is not well covered by the data and it appears that 
a high value of parallax tends to change substantially the time of 
the caustic entry inside E\. Secondly, parallax effects are partly 
degenerate with lens orbital motion, and the available data sets 
are not sufficient to disentangle these two effects. 

We therefore adopt the ESBL parameters ( s , q, uq, a, t E , to) 
found for model ESBL±P±LOM and uq > 0. Following the ar¬ 
guments summarised in the previous paragraph, other parame¬ 
ters such as parallax or orbital motion parameters cannot be used 
to constrain the lens mass and distance. 


4. Physical parameters 

4.1. Mass-distance relation from the source radius 


Combining Eq. 0 and Eq. 0 yields the following mass- 
distance relation for the lens 


M = 


P 2 K 7T re 1 


( 8 ) 


where M is the lens mass and 7r re i is a proxy to lens distance Z>l 
according to Eq. 0- In this equation, p is measured from the 
light curve modelling (Sect. 3.6 1 and 0s is derived below. 

To characterise the source, we use a DoPHOT-based data re¬ 
duction ( [Schechter et al.|1993] > of Danish I and V data to build 
an instrumental (/, (V-/)) CMD (Fig. [4]). Fitting the /-band data 
for the best-fit model parameters provides / S DK = 19.85 ± 0.01. 
We use the RCG method to measure the interstellar extinction. 
Red clump giant stars are found to be relatively sparse, which 
makes it difficult to locate precisely the mean instrumental RCG 
position, (Zrcg.dk, (V- Z)rcg,dk)- To overcome this problem, we 
cross-calibrate the Danish CMD with a CMD obtained at CTIO 
in a very similar I filter, by cross-identifying a few clump stars. 
We then measure /rcg.dk = 15.09 and (V-I) rcg.dk = 0.11. The 
values of M/.rcg = -0.12 ± 0.09 and (V - /)rcg,o = 1.06 ± 0.12 


are taken from |Nataf et ah] ( |2013| > at the coordinates of MOA- 
2007-BLG-197. Measuring the dereddened apparent magnitude 
of the red clump stars located in the Galactic bulge, |Nataf et al. | 


( |2013| > found a distance to the Galactic centre /J ( ,c = 8.20 kpc 
and an angle between the Galactic bar and the line of sight from 
the Sun </> = 40°. From these values, the distance to the RCG is 


Drcg = 


P GC sin 4> 
cos ( b ) sin (/ + <p ) 


(9) 


For the Galactic coordinates of MOA-2007-BLG-197, the RCG 
is in the far side of the bar at a distance Drcg = 8.29 kpc. We 
adopt a source at the same distance, corresponding to a distance 
modulus of p = 14.6±0.3 mag, where the error bars also account 
for the uncertainty in the position of the RCG. 

We then calculate / s ,o = Zs.dk + M/,rcg + P ~ /rcg.dk and 
obtain Is.o = 19.24 ± 0.31, which gives an absolute magnitude 
of Mi = 4.65 ± 0.31. We use a lOGyr and solar metallicity 
isochrone from Bressan e t ah] ( |2012] > to get the corresponding 
absolute magnitude My of the source from which we derive an 
intrinsic source colour of (V-Z)s.o = 0.88±0.1. This value can be 
compared to an independent estimate of the source colour based 
on the method of |Gould et al.| ( 2010a| >. MOA R^-band images 
are reduced with DoPHOT, and stars are cross-matched between 
Danish /, V and MOA R M frames (lower panel of Fig. [4]). We 
obtain 


7s,dk~7s,moa - -(0. 950±0.006)-(0.21±0.02)x(V-/) s ,dk , (10) 

where Zs.moa = 20.76 ± 0.02 from the fit of the Rm light curve 
using the best-fit parameters. Hence, Zs.dk - Zs.moa = -0.91 ± 
0.02, which yields (V - Z)s,dk = —0.18 ± 0.10. It follows that 
the estimated de-reddened source colour is (V - Z)s,o = (V — 
Z)s,dk + (V - I )rcg ~ (V - Z)rcg.dk = 0.77 ± 0.2, which is in 
good agreement with the previous estimate. In the following, we 
adopt (V - Z)s,o = 0.88 ±0.1 (most robust estimate) and Zs.o = 
19.24 ± 0.31. The source is thus a G6-K0 Main Sequence star. 

From |Kervella & Fouque|p008| brightness-colour relations, 
we estimate the angular radius of the source, 

log(fls) = 3.198-0.2/s,o±0.4895(V-/)s.o-0.0657(V-/) 2 0 , (11) 

which yields 0s = 0.54 ± 0.05 pas. With a source at Z>s = 8.3 
kpc, the physical source radius is R$ = Z>s0s = 0.96 R Q . From 
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Fig. 4. Upper panel: (/, V - I) Danish instrumental CMD of the field 
around MOA-2007-BLG-197, not corrected for interstellar extinction. 
The red point marks the position of the RCG. The magnitude of the 
microlensing target is in blue (i.e. magnified source, lens and blended 
light). The orange point corresponds to the source alone. Lower panel: 
empirical linear relation between instrumental colours (/ DK - /moa) and 
(V - I) dk (ted and green points are astrometric matches that are un¬ 
crowded in the MOA image; after sigma-clipping, only red points are 
used for the fit). 


Carroll & Ostlie (2006|), this radius is that of G8 Main Sequence 


star, which is in very good agreement with the constraint from 
the colour. 

The lens mass-distance relation obtained from Eq. 0 where 
6s = 0.54 yuas and p - 5.30 x 1(T 4 is plotted in Fig.p] 


4.2. Mass-distance relation from NaCo data 

Calibrated (J, H, K s ) NaCo magnitudes (Sect. |2.3| ) provide inde¬ 
pendently from Eq. ([8]i further lens mass-distance relations (writ¬ 
ten below for J only) through 


mj(L) = Mj + 5 log D l - 5 + Aj, 


Article number, page 8 of 16 


( 12 ) 



(J - Ks) 0 


Fig. 5. Colour-magnitude diagram {Mj, (J-K s ) o) of the lens. The curves 
are a set of isochrones from 2 Gyr (in pink) to 8 Gyr (in blue) with solar 
metallicity i jBressan et al.|2012j ). The black cross indicates the range of 
colours and magnitudes explored by the MCMC at the l-<x level. 


where m j(L) is the lens apparent reddened magnitude of the lens, 
Aj the interstellar absorption (given in Sect. |2.3| ), and Mj the 
lens absolute magnitude (which is a function of its mass M). 
The source (J, II, K s ) magnitudes are expressed as (e.g., for J ) 

mj(S) = Mj + 5 log Ds - 5 + Aj - 2.5 log(A), (13) 


where A = 1.3 is the magnification of the source at the time of 
the NaCo observations, and assuming that the absorption mostly 
occurs in the first kiloparsecs towards the bulge. The resulting 
(J, H, K s ) magnitudes are then computed from the previous equa¬ 
tions. 

The absolute magnitudes M JHK are computed using the 
isochrones by Bressan et al. ( 2012| > for both the source and lens. 
For the source we use the same isochrone as in Sect. 14.11 For 
the lens, we assume an age of 2.5 Gyr and solar metallicity. In 
fact, as seen in Fig. [ 5 ] the lens lies in a region of the isochrones 
where the age has only very small impact on the magnitudes. 
The choice of metallicity has slightly more impact, but keeps 
the variation lower than the magnitude error bars for metallicity 
changes as large as Z 0 ± 0.02. 

We then use a MCMC algorithm (Sect. |4~4{ > to match the ob¬ 
served NaCo magnitudes listed in Table [2] (Epoch 1) to the the¬ 
oretical magnitudes, using the lens mass M and distance D\ as 
fitting parameters. The observed magnitudes serve as Gaussian 
prior distributions with standard deviation equal to the magni¬ 
tudes error bars; they are shown as filled curves in Fig. [ 7 ] In the 
same figure, posterior distributions for the resulting magnitudes 
(lens, and source+lens) are shown in solid and dashed lines re¬ 
spectively. The final lens mass-distance relations are plotted in 
Fig. [6] It can be seen that while in principle three relations are 
obtained, they are very strongly correlated and thus only corre¬ 
spond to one effective independent mass-distance relation. 


4.3. Mass-distance relation from parallax or orbital motion 

In principle, when parallax is measured, a further independent 
lens mass-distance relation can be obtained by combining Eq. 0 
with Eq. 0’ 

M = , (14) 

where 7 Te is derived from the light curve modelling. Neverthe¬ 
less, as discussed in Sect. |3.6| the parallax is degenerate with 
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D l (kpc) 


Fig. 6. Lens mass-distance relations derived from the source size p (1-cr 
around the best value shown by the green shadow limited by the dashed 
lines) and NaCo ( J , H, K s ) colours constrains (respectively yellow, or¬ 
ange, red from top to bottom below the label “colours”). The blue con¬ 
tours (1 - 4 tr) represent the joint posterior probability density P(M, DO 
from the MCMC run. 



Fig. 7. Probability densities of NaCo ( J,H,K S ) magnitudes, not cor¬ 
rected for interstellar extinction. Filled curves are Bayesian Gaussian 
priors based on the NaCo measurements, which are compared to pos¬ 
terior distributions (solid lines). The dotted posterior distributions are 
those of the lens magnitudes alone (i.e., disentangled from the source; 
respectivement J, H, K s from the right to the left). 


the lens orbital motion, which prevent us from using this con¬ 
straint. We however plot on Fig. [ 6 ] the typical behaviour of the 
mass-distance relation obtained from parallax measurements, if 
7 Te had been measured. We discuss this problem further below. 


4.4. Combination of the constraints and physical parameters 


We combined the different mass-distance relations discussed in 
the previous sections (and shown in Fig. | 6 ]> in a MCMC min¬ 
imisation process to recover the lens mass M and the observer- 
lens distance Dl- We use the measurement and error bars for p, 
(J, II, K s ) and n E to build up Gaussian priors for the Bayesian 
analysis, and flat (uninformative) priors for M and Dl- The val¬ 
ues of the parameters correspond to the solution uq > 0 , as 
discussed in Sect. |3.6| We used the convergence criterion of 


Geweke (1992) to stop the MCMC and compute the posterior 


probability densities. 

As seen in Fig. [ 6 ] in principle this problem is over¬ 
constrained, with three relations for only two fitting parame¬ 
ters. This enabled us to check the consistency of the different 


measurements, and investigate further the degeneracy between 
parallax and lens orbital motion, which prevents n e from being 
used as a constraint. We first noticed a clear discrepancy between 
the value of tte derived from the fit (the prior) and the poste¬ 
rior value. This pointed out a strong tension between the three 
constraints from parallax, source size and infrared colours. We 
then successively removed the parallax, source size or colours 
constraints from the MCMC runs. Without the source size con¬ 
straint, a model with parallax values as large as 7 Te ~ 1.1 was 
found. This would require a much closer lens (Dl ~ 0.6 kpc), a 
longer Einstein radius (6 E ~ 1 .6 mas) and consequently a much 
smaller source size (p ~ 3.3 x 10 4 ). This value of p is many 
sigma away from that measured from the light curve modelling. 
Such a difference is very unlikely given the very strong con¬ 
straint on p obtained from the caustic exit modelling. This there¬ 
fore confirmed previous concerns that in this case, the degen¬ 
eracy between parallax and lens orbital motion does not allow 
us to use the parallax mass-distance relation. Consequently, the 
parallax constraints are removed from the MCMC to derived the 
final values for M and Dl, thus constrained by the source size 
and the colours. 

The final values for M and Dl are given in Table [4] with the 
correlation plot shown in Fig. [ 6 ] We also compute the mass of 
the lens companion, m BD = qM, and find ;n BD = 41 ± 2Mj. 
The host mass is 0.82 ± 0.04 Mq. The companion is thus a brown 
dwarf orbiting a solar-type star. 


Physical parameter [unit] 

Value 

Lens mass, M [M 0 ] 

0.86 ± 0.04 

Host mass [M Q ] 

0.82 ± 0.04 

Companion mass, ;hbd [Mj] 

41+2 

Projected orbit, a ± [AU] 

4.29 + 0.10 

Observer-lens distance, Dl [kpc] 

4.17 + 0.30 

D s [kpc] 

8.3+ 1.0 

6 S [//as] 

0.54 ± 0.05 

0 E [mas] 

0.91+0.04 

v ± [km.s -1 ] 

80 + 2 

7Te (calculated) [mas] 

0.13 + 0.01 

y (new fit) [year -1 ] 

0.26 ± 0.07 

P 

0.08 ± 0.07 


Table 4. Physical parameters of MOA-2007-BLG-197L and its com¬ 
panion. 


We finally perform additional consistency checks on the 
overall parameters of MOA-2007-BLG-197L. First, we compute 
the projected lens-source relative velocity 


v± = 




(15) 


and find v ± - 80km.s , which is in very good agreement with 
probability densities predicted by |Dominik] ( ]2006[ > for the de¬ 
rived values of t E and Dl- Moreover and as a supplementary 
check, we computed the probability distribution of 7 Te from the 
distributions of M and Dl using Eq. ( fl4j >, and derived its maxi¬ 
mum a posteriori (MAP) value. We then ran again the light curve 
modelling with this value of n E kept fixed. For the overall consis¬ 
tency, all the parameters are fixed as well, except s and a (values 
are given in Table [3] model with uq > 0) which yields the value 
of the orbital motion parameter y (see Table |4|. From this, we 
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compute 


P = 


Fkin.l 


-'pOt.-L 


2AU 2 


3,,2 


7T E s-'y 


c Oe (^e + ks/9e) 


(16) 


the ratio of the apparent kinetic to potential energy for the bi¬ 
nary lens orbit projected onto the plane of the sky (An et al. 
|2002) >. Gravitationally bound systems should exhibit 0 < /? < 1 
(for face-on circular orbits, p = 0.5). Using = AU/Ds, we 
find p =* 0.08, which is a value consistent with a bound sys¬ 
tem. In contrast, high parallax values of e.g., Ue = 1.1 lead to 
much lower values of p ^ 4.5 x 10~ 3 . As argued by Batista et al. 


( j201l| , a scenario with /? <s 1 would require highly improbable 
orbital parameters or projections, such as a very close-in com¬ 
panion seen on a nearly edge-on circular orbit, which is excluded 
here. 


5. MOA-2007-BLG-197Lb in the brown dwarf 
landscape 

While around 1900 exoplanets have been detected so fai[^] less 
than a hundred brown dwarfs orbiting stars are known today. The 
lack of brown dwarf companions to solar-type stars was noticed 
early after the first exoplanets detections, and referred to as a 
brown dwarf desert”. Using the radial velocity method, |Marcy| 


& Butler (2000) noted that at tight orbital separations (a < 3AU), 


companions more massive than 8 Mj represented only a very low 
fraction of the detected objects. This was not expected from the 
relatively high frequency (~ 13%) of binary stars found earlier 
by Duquennoy & Mayor| ( fl991) l, also at close separation. 

This lack of brown dwarf companions compared to plane¬ 
tary companions has since been confirmed by all detection meth¬ 
ods. For instance, Lafreniere et al. (2007) found a frequency 


of 1.9*H% for 13-75 Mj companions located in the range 
25 - 250 AU from the Gemini Deep Planet Survey. Similarly, 
Metchev & Hillenbrand ( 2009| derived a frequency of3.2±3.1% 


for brown dwarfs orbiting young solar-type stars in the range 


28 - 1590 AU using adaptive optics direct imaging. McCarthy & 


Zuckerman (2004) also found a higher proportion of stellar com¬ 


pared to substellar companions to stars, and derived a frequency 
ratio of ~ 3 - 10 at wide separations. |Luhman et al. ( }2007[ > found 
that this ratio is comparable to the relative abundance of stars 
to substellar objects when they are found single, either in star 
forming regions or in the solar neighbourhood. 

The upper panel of Fig. [8] is a mass vs. period diagram 
which diplays the known brown dwarf companions of solar- 
type stars detected through microlensing, radial velocimetry, 
transit and direct imaging. For radial velocimetry and tran¬ 
sit, we included objects from Ma & Ge ( 2014[ ) catalogue with 
(minimum) masses in the range 13 - 74 Mj and orbital peri¬ 
ods below 2 x 10 4 d. Furthermore, we have excluded objects 
with mass uncertainties above 25 Mj, and deleted (false posi¬ 
tive) TYC 1240-945-1. Microlensing brown dwarfs are shown as 
green points, they are: MOA-201 l-BLG-149 and OGLE-2011- 
BLG-0172/MOA-201 l-BLG-104 (JShin et al.||2012bj), MOA- 
2010-BLG-073 dStreet et al.| |2013| i. OGLET’013-BLG-0578 
( Park et al.|2015[ ) as well as MOA-2007-BLG-197 (this paper). 
For consistency with Doppler and transit data, we did not in¬ 
clude MOA-2009-BLG-411 ( jBachelet et aT|f20T2] ) which has a 
large uncertainty in the mass. In the case of microlensing objects 
for which only projected separations a ± were measured (all ob¬ 
jects in this case), we estimate the (physical) semi-major axis a 
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Fig. 8. The upper panel shows a mass-period diagram with brown 
dwarfs companions detected through radial velocity, transit and direct 
imaging (blue filled circles for measured masses, and open circles for 
minimum masses; list adapted from |Ma & Ge|20l4j >, and microlensing 
(green diamonds). For reference, exoplanets detected so far are also dis¬ 
played (small black circles, http://exoplanets.org). The red dot¬ 
ted line indicates the global radial velocity completeness limit, while 
the red dashed line marks the region above which data are included to 
perform the non-parametric, two-dimensional probability density distri¬ 
bution shown in the lower panel. 


as the median of the probability distribution of Gould & Loeb 
([1992), which we find to be equal to a = (2/ V3 )a± « 1.15aj_. 
Kepler’s third law is then used to yield the estimated period. The 
error bar on measured a ± is propagated to a following the previ¬ 
ous equation, and the error bar on the period is obtained through 
a MCMC run including uncertainties on a and primary mass M 
(Gaussian distributions are assumed). 

As seen in Fig. [8] the distribution of brown dwarfs is not uni¬ 
form, in particular it exhibits an increasing frequency of objects 
with increasing orbital period. Furthermore, |Ma & Ge ( 2014| > 
argue that the driest part of the brown dwarf desert seems to 
be confined within a region at close separation, namely within 
P < lOOd and for masses within 30 - 55 Mj for Sun-like host 
stars. They subsequently split the objects into two distinct popu¬ 
lations, whether their mass is higher or lower than 42.5Mj: the 
more massive brown dwarfs would mainly be the outcome of 
gravitational fragmentation and collapse of a molecular cloud 
(star-like formation scenario, corroborated by an eccentricity 
distribution similar to binary stars), while brown dwarfs below 
42.5 Mj would mainly result from gravitational instability within 
the proto-planetary disk. But other authors such as[Guillot et al. 


( 2012| >, however, argue that the depletion of objects at tight or¬ 
bits may as well be explained by a loss of an initial population 
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of close-in brown dwarfs due to tidal interactions with their host 
stars. In this scenario, close-in, massive objects lose angular mo¬ 
mentum due to the slower rotation of the star relative to the plan¬ 
ets orbital motion, spiral in and fall into the star. This effect is 
predicted to peak for G-dwarf primaries and should not be im¬ 
portant for earlier-type stars, which is supported by the detection 
of close-in brown dwarfs around F-type stars. 

Different mechanisms may overlap to shape the brown 
dwarfs landscape, and lead to a more complex structure than 
previously discussed. To analyse it further, we performed a non- 
parametric, two-dimensional probability density distribution fit 
to the data, only for objects above the red dashed line (i.e., above 
the radial velocity completeness limit, based on Mayor et al. 
( 2011| , and marked by the red dotted line in Fig. [8}. The prob¬ 
ability density estimation is based on 
we used a Gaussian kernel, and follow 


Scott (1992 1 ; in practice 
Silverman| ( |1986 1 rule to 


estimate the bandwidth parameter. The resulting probability den¬ 
sity is shown in the lower panel of Fig. [8] We used two different 
methods to compute the bandwidth parameter, Silverman (|1986|l 


Scott ((1992 1 , and found that the resulting probability density 


is relatively stable, even if a few data points are removed from 
sparsely sampled regions. 

The first striking feature of the density profile (lower panel 
of Fig. [8j is a region of depletion of objects at intermediate 
masses {M ~ 30 - 60 Mj) and short orbital periods (P < 30 d). 
This matches the region referred to as the driest part of the 
brown dwarf desert according to |Ma & Ge| ( |2014[ ). Second, an 
accumulation of objects can be seen around P ~ 500 d and 
M ~ 20 Mj (following the apparent trend observed for giant 
planets, black data points in Fig. [5). Third, we find another de¬ 
pletion of objects at long orbital periods (P > 500 d) and high 
masses (M > 50 Mj). Since the brown dwarf sample is drawn 
from various surveys (mainly from radial velocities), the inter¬ 
pretation of these features should be taken with caution. Al¬ 
though all brown dwarfs in Fig.[8]are chosen above the radial ve¬ 
locity completeness limit from jMayor et al.| ( j2011) , it is not guar¬ 
anteed that these objects are not affected by observational biases. 
In particular, a degeneracy between mass and period is increas¬ 
ingly affecting very long-period brown dwarfs, and the different 
surveys should be corrected from their sample size. Hence, while 
a thorough analysis of all these factors is necessary to assess the 
exact shape of the brown dwarfs distribution, the gross features 
emerging from the lower panel of Fig.[8]may well be real. If we 
split the density profile into two regions of masses above and be¬ 
low 42.5 Mj as discussed by Ma & Ge] ( |2014[ ), the distribution of 
high-mass objects can be seen as shifted towards shorts orbital 
periods (P ~ 30 d), while less massive objects appear to pile up 
at longer orbital periods (P ~ 500 d). This supports the claim 
of Ma & Ge| ( 2014j l that massive objects would accumulate at 
short periods as a result of gravitational collapse of a molecular 
cloud, contrary to less-massive objects built up by gravitational 
instability in the disk which would accumulate at longer periods. 
Nevertheless, the depletion at intermediate masses and short pe¬ 
riods would as well be affected by a star-engulfing mechanism 
advocated by Guillot et al. ( |2012| ), in particular because almost 
all brown dwarfs included here orbit Sun-like stars, for which 
this effect peaks. 

The distribution of brown dwarfs as a function of mass and 
period still remains uncertain because of a lack of detections. 
It is thus difficult to distinguish between the different mecha¬ 
nisms which shape the brown dwarf desert. Moreover, these de¬ 
tections are mostly objects orbiting Sun-like stars, which makes 
it difficult to study their distribution as a function of host star 
mass. While the MOA-2007-BLG-197L companion orbits a G-K 


dwarf host, microlensing hosts are most frequently M-K dwarfs. 
Hence, future microlensing surveys will provide unique infor¬ 
mation on the brown dwarf distribution for low-mass hosts, thus 
offering a complementary mass bin to other methods. 


6. Summary and prospects 


We have presented MOA-2007-BLG-197Lb, the first brown 
dwarf discovered around a Sun-like star through gravitational 
microlensing. The system is located at 4.2 + 0.3 kpc from the 
Earth, the brown dwarf companion has a mass of 41 + 2Mj and 
was observed at a projected distance of 4.3 + 0.1 AU from a K 
dwarf host star. We have performed a non-parametric probabil¬ 
ity density distribution fit to the population of brown dwarf com¬ 
panions detected so far. The resulting density profile has a com¬ 
plex structure, leading to a “brown dwarf landscape” most likely 
shaped by different (and perhaps competitive) formation or de¬ 
struction mechanisms. While it seems difficult with the current 
data set to distinguish observationally which are the dominant 
mechanisms, an answer to this question appears to within the 
reach of further observations in the near to mid-term future. 

Gravitational microlensing is an exceptional tool to detect 
brown dwarfs as free-floating objects, companions to stars or as 
brown dwarfs binaries. It has a unique sensitivity to detect brown 
dwarfs companions to stars of any type, in particular at long or¬ 
bital periods. Recent advances in using networks of robotic tele¬ 
scopes (e.g., [Tsapras et al7| |2009| will provide in a near future 
an order of magnitudes more brown dwarfs detections through 
microlensing. Future microlensing space missions ( |Penny et al.| 
2013||Yee et al.|2014| Beaulieu et al.|2010|> also carry important 


promises for providing unique information on the populations of 
brown dwarfs in their different configurations. 
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Fig. 9. Correlations between the parameters derived from ESBL’s model discussed in Sect. |3.l| The red point refers to the best-fitting model 
obtained. 
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Fig. 10. Idem for ESBL+P's model. 
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Fig. 11. Idem for ESBL+LOM’s model. 
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Fig. 12. Idem for ESBL+P+LOM"s model ( u 0 > 0) with p' = 5.3 x 10 4 . 


Article number, page 16 of 16 


CDF 
















































